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Abstract 

The Kadomtsev-Petviashvili (KP) equation describes weakly dispersive and small amplitude waves propagating in 
a quasi-two dimensional situation. Recently a large variety of exact soliton solutions of the KP equation has been 
found and classified. Those soliton solutions are localized along certain lines in a two-dimensional plane and decay 
exponentially everywhere else, and they are called line-soliton solutions in this paper. The classification is based on 
the far- field patterns of the solutions which consist of a finite number of line-solitons. In this paper, we study the 
initial value problem of the KP equation with V- and X-shape initial waves consisting of two distinct line-solitons by 
means of the direct numerical simulation. We then show that the solution converges asymptotically to some of those 
exact soliton solutions. The convergence is in a locally defined L^-sense. The initial wave patterns considered in this 
paper are related to the rogue waves generated by nonlinear wave interactions in shallow water wave problem. 

Keywords: Kadomtsev-Petviashvili equation, soliton solutions, chord diagrams, pseudo-spectral method, window 
technique 



1. Introduction 

The KdV equation may be obtained in the leading order approximation of an asymptotic perturbation theory for 
one-dimensional nonlinear waves under the assumptions of weak nonlinearity (small amplitude) and weak dispersion 
(long waves). The initial value problem of the KdV equation has been extensively studied by means of the method of 
inverse scattering transform (1ST). It is then well-known that a general initial data decaying rapidly for large spatial 
variable evolves into a sum of individual solitons and some weak dispersive wave trains separated away from solitons 
(see for examples, im[T3l[T5ll20ll). 

In 1970, Kadomtsev and Petviashvili I?! proposed a two-dimensional dispersive wave equation to study the stabil- 
ity of one soliton solution of the KdV equation under the influence of weak transversal perturbations. This equation 
is now referred to as the KP equation, and considered to be a prototype of the integrable nonlinear dispersive wave 
equations in two dimensions. The KP equation can be also represented in the Lax form, that is, there exists a pair of 
linear equations associated with the eigenvalue problem and the evolution of the eigenfunctions. However, unlike the 
case of the KdV equation, the method of 1ST based on the pair of linear equations does not seem to provide a practical 
method for the initial value problem with non-periodic waves considered in this paper. At the present time, there is no 
feasible analytic method to solve the initial value problem of the KP equation with initial waves having line-solitons 
in the far field. 

In this paper, we study this type of the initial value problem of the KP equation by means of the direct numerical 
simulation. In particular, we consider the following two cases of the initial waves: In the first case, the initial wave 
consists of two semi-infinite line-solitons forming a V-shape pattern, and in the second case, the initial wave is given 
by a linear combination of two infinite line-solitons forming X-shape. Those initial waves have been considered in 
the study of the generation of large amplitude waves in shallow water lfT6l [TSl fuyi . The main result of this paper is 
to show that the solutions of the initial value problem with those initial waves asymptotically converge to some of the 
exact soliton solutions found in ||2l[3l. This implies a separation of the (exact) soliton solution from the dispersive 
radiations in the similar manner as in the KdV case. 
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The paper is organized as follows: In Section |2] we provide a brief summary of the soliton solutions of the KP 
equation and the classification theorem obtained in (T. T] for those soliton solutions as a background necessary for 
the present study. In particular, we introduce the parametrization of each soliton solution with a chord diagram which 
represents a derangement of the permutation group, i.e. permutation without fixed point. In Section [3] we present 
several exact soliton solutions, and describe some properties of those solutions. Each of those soliton solutions has 
numbers of line-solitons in a far field on the two-dimensional plane, say in y » 0, and numbers of line-solitons in 
the far field of the opposite side, i.e. y <sc 0. This type of soliton solution is referred to as an (A^_, A^+)-soliton solution. 
Here we consider those solitons with +A^+ < 4 and < 3. In Section]?] we describe the numerical scheme used in 
this paper, which is based on the pseudo-spectral method combined with the window technique 1.17. ,19J . The window 
technique is especially used to compute our non-periodic problem which is essentially an infinite domain problem. 
Finally, in Section|5] we present the numerical results of the initial value problems with V- and X-shape initial waves, 
and show that the solutions asymptotically converge to some of those exact solutions discussed in Section 13] The 

convergence is in the sense of locally defined L2-sense with the usual norm, i.e. •= (/Jq \f(x,y)\'^dxdyy , 

where D c is a compact set which covers the main structure describing the (resonant) interactions in the solution. 
We also propose a method to identify an exact solution for a given initial wave with V- or X-shape pattern based on 
the chord diagrams introduced in the classification theory. 



2. Background 



Here we give a brief summary of the recent result of the classification theorem for soliton solutions of the KP 
equation (see (|9] |2l El for the details). In particular, each soliton solution is then parametrized by a chord diagram 
which represents a unique element of the permutation group. We use this parametrization throughout the paper. 

2.1. The KP equation 

The KP equation is a two-dimensional nonlinear dispersive wave equation given by 



dx (4d,u + 6udxU + dluj + 3dju = 0, 



where d".u 



d"u 
dx" 



etc. Let us express the solution in the form. 



u{x, y,f) -2 d\ In t(x, y, t) . 



(2.1) 



(2.2) 



where the function t is called the tau function, which plays a central role in the KP theory. In this paper, we consider 
a class of the solutions, where each solution can be expressed by T{x,y,i) in the Wronskian determinant form r - 
Wr(/i,...,.Mi.e. 

/i /; ■■■ /i 

h /2 ■■■ f. 



with/y' := d'xfn forn = 1, 
linear equations. 



T(x,y, f) = 
, A^. Here the functions {/i , . . 



i'N-V) 



JN 



(2.3) 



, /at) form a set of linearly independent solutions of the 



dyfn = djn, d,fn = -5^./„. 



(The fact that \2.2\ with \2.3\ gives a solution of the KP equation is well-known and the proof can be found in several 

(2.4) 



places, e.g. see Il6l|3|.) The solution of those equations can be expressed in the Fourier transform, 

fn{x,yj)^ J e'''^'''>-''''p„(k)dk, « = 1,2,...,M 



2 



with an appropriate contour C in C and the measure p„(k) dk. In particular, we consider a finite dimensional solution 
with Pn{k) dk - Yjm=\ ^n,m^ik - k,„) dk with fl„_,„ e M, i.e. 

M 

f„(x, y, f) = ^ a„^,„E,„(x, y, t) with £,„ = exp(^;„x + klj - k^t). 



m=l 



kM) and the N x M matrix A :- (a„,„) of 

(2.5) 



Thus this type of solution is characterized by the parameters {ki,k2, 
rank(A) - N, that is, we have 

(/i,/2,...,/w) = (£i,£2,...,£m)A^ 

Note that {£^,£2, ■ ■ ■ ,Em} gives a basis of M*' and {/i,/2, . . . ,/a?) spans an A^-dimensional subspace of M*^. This 
means that the A-matrix can be identified as a point on the real Grassmann manifold Gr(A^, M) (see ||9] [S]). More 
precisely, let Mnxm(^) be the set of all x M matrices of rank A^. Then Gr(A^, M) can be expressed as 

Gr(A?,M) = GW(M)\M^vxm(M), 

where GLAr(IR) is the general linear group of rank A^. This is saying that other basis ig\, . . . ,gN) = if\, - ■ ■ ,fN)H for 
any H e GL^(K) spans the same subspace. Notice here that the freedom in the A-matrix with GL^(K) can be fixed by 
expressing A in the reduced row echelon form (RREF). We then assume throughout this paper that the A-matrix is in 
the RREF, and show that the A-matrix plays a crucial role in our discussion on the asymptotic behavior of the initial 
value problem. 

Now using the Binet-Cauchy Lemma for the determinant, the r-function of (|2.2|l can be expressed in the form. 



T - 



El 
kiEi 

k^'Ei 



E2 

kl^E2 ■■ 



kmEm 



(an ail 
a\2 022 



VfllM a2M 
■ , jN)E(Ji, j2, ■ ■ ■,Ef^), 



ONI \ 



(2.6) 



where ^(;i, 
given by 



l<jl<j2<-<jN<M 

, Jn) is the N X N minor of the A-matrix with columns marked by (ji,. . . , j^), and E{ji , . . ., Jn) is 



E(ji 



, M = Wr(£,, ,...,£,,) = Y\{kj^, - kj,)Ej, ■ ■ ■ Ej^ 



l<m 



From the formula (2.6 1, one can see that for a given A-matrix, 



• if a column of A has only zero elements, then the exponential term E„, with m being the column index never 
appear in the r-function, and 

• if a row of A has only the pivot as non-zero element, then E„ with n being the row index can be factored out 
from the r-function, i.e. E„ has no contribution to the solution. 

We then say that an A-matrix with no such cases is irreducible, because the r-function with reducible matrix can be 
obtained by a matrix with smaller size. 

We are also interested in non-singular solutions. Since the solution is given by M = 25^(lnr), the non-singular 
solutions are obtained by imposing the non-negativity condition on the minors, 

^(;i , ;2, • • • , Jn) > 0, for all 1 < ;i < 72 < ■ ■ ■ < ;w < M. (2.7) 

This condition is not only sufficient but also necessary for the non-singularity of the solution. We call a matrix having 
the condition (2.7 1 totally non-negative matrix. 
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2.2. Line-soliton solution and the notations 

Let us here present the simplest solution, called one line-soliton solution, and introduce several notations to de- 
scribe the solution. A line-soliton solution is obtained by a r-function with two exponential terms, i.e. the case N - 1 
and M - lin (|2.5|i: With the A-matrix of the form A = (1 a), we have 



T = Ei+aE2 = 2 VaeJ***'^*-^ cosh ^(6i - 6*2 - In a) , 

The parameter a in the A-matrix must be a > for a non-singular solution (i.e. totally non-negative A-matrix), and 
it determines the location of the soliton solution. Since a = leads to a trivial solution, we consider only a > (i.e. 
irreducible A-matrix) . Then the solution u = 25 ^(In t) gives 

1 9 9 1 

u - -(k\ - ^2) sech -{6\ -62- In a). 

Thus the solution is located along the line 9\ - 02 - In a. We here emphasize that the line-soliton appears at the 
boundary of two regions, in each of which either E\ or E2 becomes the dominant exponential term, and because of 
this we also call this soliton [l,2]-soliton solution (or soliton of [l,2]-type). In general, the line-soliton solution of 
[/, 7] -type with / < j has the following structure (sometimes we consider only locally), 

u = A[,-,,.] sech^ i (K[,,,j ■ X - ^[yjf + 0[,;,]) (2.8) 
with some constant 0[,j]. The amplitude A[,j], the wave-vector K[,j] and the frequency are defined by 

K[,,] = [kj - ki, k] - kf) = (kj - fc,.) (1 , k^ + kj) , 
= - k] = (kj - ki)(kf + kikj + k]). 

The direction of the wave- vector K[,j] = (K^. .^K^. is measured in the counterclockwise from the jc-axis, and it is 
given by 

''•^^ -tan>P[,j] ^ki + kj. 



Kf. , 



Notice that ^[ij] gives the angle between the line K[, j] ■ x = const and the y-axis (See Figure [2?T] l. Then a line-soliton 
(2.8 I can be written in the form with three parameters A *!'[,■ and x'-^j jy 



u = A[y] sech^ A^^l^ [x+y tan 'Piij] - C[,j]f - 4,/]) . (2-9) 

with C[ij] = kf -i- kikj + kj - jAyj] + | tan^*I'[,j]. For the parameter x^.^.j giving the location of the line-soliton, we 
also use the notation. 



with in (2.8 1 which is determined by the A-matrix and the A: -parameters. For multi-soliton solutions, one can 
only define the location of each [/, y]-soliton using the asymptotic position in the xy-plane either x » or x <K 
(we mainly consider the cases where the solitons are not parallel to the y-axis), and we use the notation x'^^j-^ (or x'^^j-^) 
which describes the x-intercept of the line determined by the wave crest of [/, j] -soliton in the region x » (or x <K (3) 



at f = 0. In Figure 2. 1 we illustrate an example of one line-soliton solution. In the right panel of this figure, we 
show a chord diagram which represents this soliton solution. Here the chord diagram indicates the permutation of the 
dominant exponential terms and Ej in the r-function, that is, with the ordering A:, < kj, E, dominates in x <K 0, 
while Ej dominates in x » 0. This representation of the line-solitons in terms of the chord diagrams is the key concept 



throughout the present paper (See section 2.3 below for the precise definition of the chord diagrams.) 
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Figure 2.1: One line-soliton solution of [i,j]-type and the corresponding chord diagram. The amplitude and the angle ^[ij] are given by 
A[ij] = 1(4/ - kj)^ and tan'P[, j] = kj + kj. The upper oriented chord represents the part of [i, y]-soliton for y » and the lower one for )' <K 0. 
Here /, j represent not just kj, kj, but also the value of those parameters, i.e. the horizontal line is the coordinate for the i:-parameters, and in this 
example, ki = -0.5 and k2 = I. 



For each soliton solution of (2.9i, the wave vector K[, j] and the frequency Q[,j] satisfy the soliton-dispersion 
relation, i.e., 

A^[,j,Kl..^^{Kl,J + 3{KlJ. (2.11) 
The soliton velocity V[; ,] defined by K;,- ,] ■ V[, ,] = ,] is given by 



Qr- 1 k- + kikj + kj 



|Kr,- 



1 + {ki + kjY 



Note in particular that Cyj] - kf + kikj + kj > 0, and this implies that the x-component of the velocity is always 
positive, that is, any soliton propagates in the positive jc-direction. On the other hand, one should note that any small 
perturbation propagates in the negative x-direction, that is, the .^-component of the group velocity is always negative. 
This can be seen from the dispersion relation of the KP equation for a linear wave (p - exp(/k ■ x - ia>t) with the 
wave-vector k = {k^, ky) and the frequency cj, 

1,3 3^? 

from which the group velocity of the wave is given by 





d(jj \ 


' 3 


( ]P-\ 
k^ + - 


2 k," 


[dkx 


dk,]^ 
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This is similar to the case of the KdV equation, and we expect that asymptotically soliton separates from small 
radiations. Physically this implies that soliton is a supersonic wave due to its nonlinearity (recall that the velocity 



of shallow water wave is proportional to the square root of the water depth, and the KP equation in the form (2.1 
describes the waves in the moving frame with the phase velocity in the x-direction). 

Stability of one-soliton solution was shown in the original paper by Kadomtsev and Petviashvii |7|, and this may 
be stated as follows; For any e > and any r > 0, there exists 5 > so that if the initial wave u{x,y,Q) :- u^{x,y) 
satisfies 

- «exactllL2(Dl') < ^. 

for some exact soliton solution, ul,^^^^ := Aq sech^ sj^i^ + y *J'o - Cot - xq) with appropriate constants Aq, and 
xo (recall Co = 5A0 + | tan^ *Po), then the stability implies that the solution u'{x,y) :- u(x,y, t) satisfies 

- «exactllL2(Z5;.) < for t 00, 

where D' is a circular disc with radius r moving with the soliton, i.e. 

d; = {(x,y) eR^:(x- xo(t)f + (y - yo(t)f < r\ 



with (xo(t),yo(t)) at any point on the soliton, xo(f) + y(}{t)tan^'o = Cot + xq. Here ||/||z,2(£,) is the usual L^-norm of 
f(x,y) over a compact domain D c M^, i.e. 



\\LHD)-^l^j£\fix,yfdxdy 



This stabiHty impHes a separation of the soliton from the dispersive radiations (non-soliton parts) as in the case of the 
KdV soliton. We would like to prove the similar statement for more general initial waves. However there are several 
difficulties for two-dimensional stability problem in general. In this paper, we will give a numerical study for some 
special cases where the initial waves consist of two semi-infinite line solitons with V- or X-shape. 

Finally we remark that a line-soliton having the angle *P ^ | has an infinite speed, and of course it is beyond the 
assumption of the quasi-two dimensionality. We also emphasize that the structure of the solution for y » can be 
different from that in y <K 0, that is, the set of asymptotic solitons in y » can be different from that in y <k 0. This 
difference is a consequence of the resonant interactions among solitons as we can see throughout the paper 

2.3. Classification Theorems 

Now we present the main theorems obtained in fT, ^3] for the classification of soliton solutions generated by the 
T-functions with irreducible and totally non-negative A-matrices: 

Theorem 2.1. Let {ei, . . . ,ei^} be the pivot indices, and let [gi,. . . ,gM-N] be non-pivot indices for an N x M irre- 
ducible and totally non-negative A-matrix. Then the soliton solution generated by the T-function with the A-matrix 
has the following asymptotic structure: 

(a) For y » 0, there are N line-solitons of [e„, j„]-typefor n — I, . . . ,N. 

(b) Fory <K 0, there are M — N line-solitons of[i„„g„]-typefor m — I, . . . , M — N. 
Here j„(> e,,) and im(< gm) ore determined uniquely from the A-matrix. 



The unique index pairings [e„,j„] and [im,gm] in Theorem 2.1 have a combinatorial interpretation. Let us define 
the pairing map tt : {1,2, ... , M] — > {1,2,..., M] such that 

{ 7T(e„) = j,„ n = 1,...,A^, 

\ (2.12) 
[ Mgm) = im, m= 1,...,M -N, 

where e„ and g„, are respectively the pivot and non-pivot indices of the A-matrix. Then we have: 

Theorem 2.2. The pairing map n is a bijection. That is, n e S m, where S m is the group of permutation for the index 
set [1,2,..., M}, i.e. 

lei ■■■ Cn gl ■■■ gM-N\ 

\ji ■ ■ ■ Jn h ■ ■ ■ Im-n I 



n = 



Note in particular that the corresponding n is the derangement, i.e. n has no fixed point. 



Theorem 2.2 shows that the pairing map tt for an (M - A^, A^)-soliton solution has excedances, i.e. n{i) > i for 
i - I,. . . ,N, and the excedance set is the set of pivot indices of the A-matrix. We represent each soliton solution with 
the chord diagram defined as follows: 

(i) There are M marked points on a line, each of the point corresponds to the fc-parameter. 

(ii) On the upper side of the line, there are chords (pairings), each of them connects two points on the line 
representing 7r(e„) - j„ for n - I, . . . ,N, i.e. the excedance j,, > e„. 

(iii) On the lower side of the line, there m'e M - N chords representing n{gm) = !,„ for m = 1, . . . , M - N, i.e. the 
deficiency /,„ < g^- 
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In Figure 



2.2 



we illustrate an example of the chord diagram which represents the derangement, 



or simply tt = (46523817). 

The diagram then shows that the set of excedances is {1,2,3,6), and the corresponding soli ton solution consists of 
the asymptotic line-solutions of [1,4]-, [2,6]-, [3,5]- and [6,8]-types in y » and of [1,7]-, [2,4]-, [3,5]- and 
[7,8]-types iny <K 0. 




Figure 2.2: Example of the chord diagram. This shows n = (46523817). Each chord joining /-th and y'-th marked points represents the line-soliton 
of [(', ;]-type. 




1 2 3 4 5 6 7 8\ 
4 6 5 2 3 8 1 7/ 



3. Exact solutions 

Here we present several exact solutions generated by smaller size matrices with N < 3 and M < 4. Those solutions 
give a fundamental structure of general solutions, and we will show that some of those solutions appear naturally as 
asymptotic solutions of the KP equation for certain classes of initial waves related to the rogue wave generation 
lfT6l [T8l [191 l5l . The detailed discussions and the formulae given in this section can be found in [31 • 

3.1. Y-shape solitons: Resonant solutions 

We first discuss the resonant interaction among line-solitons, which is the most important feature of the KP equa- 
tion (see e.g. [TTl T?" "81). To describe resonant solutions, let us consider the r-function with M = 3, that is, the 



T-function has three exponential terms {/ii , £2, £3 ). In terms of the r-function in the form ( 2.3 1, those are given by the 
cases (A? = 1, M = 3) and (A^ = 2, M = 3). 

Let us first study the case with N - 1 and M = 3, where the A-matrix is given by 



A = (1 a b). 



The parameters a, b in the matrix are positive constants, and the positivity implies the irreducibility and the regularity 
of the solution. The r-function is simply given by 

T - E\ + aEo + bE^. 

(Note that if one of the parameters is zero, then the r-function consists only two exponential terms and it gives one 
line-soliton solution, i.e. reducible case.) With the ordering ki < k2 < kj,, the solution u - 25^-(lnr) consists of 
[1, 3]-soliton for y » and [1,2]- and [2, 3]-solitons for y «; 0. Taking the balance between two exponential terms 
in the r-function, one can see that those line-solitons of [1, 3]- and [2, 3]-types are localized along the lines given in 



(2.9i, i.e. X + ytan^'[ij] - C[ij]t - for - [1,3] and [2,3] where the locations xj'.^.j are determined by the 
A^natrix (see ( |2.10| l), 

Xn \nb, x%^,- In-. (3.1) 

['■^1 |3,1| |3,2| a ^ 

where |/, j| := kj - kj. The shape of solution generated by r = £1 -1- aE2 + bE^ with a - b - I (i.e. at f = three 



line-solitons meet at the origin) is illustrated via the contour plot in the first row of Figure 3.1 This solution represents 
a resonant solution of three line-solitons, and the resonant condition is given by 

K[i,3] - K[i,2] + K[23], ^[1,3] = ^[1.2] + ^[2,3]> 
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T=0 




which are trivially satisfied with K[, j] = (kj - ki,kj - k^) and Q[, j] - - k\ The chord diagram corresponding to 
this soUton is shown in the right panel of the first raw in Figure [TT] and it represents the permutation n - (312). 
Let us now consider the case with N - 2 and M - 3: We take the A-matrix in the form, 

1 -b 
1 a 

where a and b are positive constants. Then the r-function is given by 

T = £(1,2) + fl£'(l, 3) + bE(2,3), 

with E{i, j) = (kj - ki)EiEj for i < j. In this case we have [1,2]- and [2, 3]-solitons for y » and [1, 3]-soliton for 
y ^0, and this solution can be labeled by tt = (231). Those line-solitons of [1,2]- and [1, 3]-types are localized along 
the lines, x -H y tan jj - C[,j]f - Jf|'.^.j with 

1 , (\i,2\b\ 1 , /I3,2| 

= - ^ 1" ( ^ « ) ' ^[>3, — ^ In ( I - (3-2) 



where |/, j\ :- kj - kj. In the lower figures of Figure 3.1 we illustrate the solution in this case. Notice that this figure 
can be obtained from (2, l)-soliton in the upper figure by changing (x,y) — » {-x, -y). Here the parameters a and b in 
the A-matrix are chosen, so that jc^ = 3] = 0, that is, all of those soliton solutions meet at the origin at f = 0. 

3.2. N - 1,3 and M = 4 cases 

Let us first discuss the case with N - 1 and M = 4, that is, the A-matrix is given by 

A - (I a b c), 

where a, b and c are positive constants. The r-function is simply written in the form 

T = El + aE2 + bEj, + CE4. 
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Figure 3.2: Examples of (3,1)- and (1, 3)-soliton solutions. The solution in the upper figures are generated by the 1x4 A-matrix with 
(ki , k2, ki,k4) = (- V3, -1,1, V3). The symmetry in the <:-parameters implies that [1,4]- and [2, 3]-solitons are parallel to the 7-axis. The lower 
figures are generated by the 3x4 A-matrix with (ki,k2, kj, k4) = ~ 2^/2' 



In this case, we have (3, l)-soliton solution consisting of one line-soliton of [l,4]-type for y » and three hne- 
solitons of [1, 2]-, [2, 3]- and [3,4]-types for y <^ Q. This solution is labeled by tt = (4 1 2 3). The upper figures in 



Figure 3.2 shows the time-evolution of the solution of this type. We seta = b = c-lof the A-matrix, so that all four 
line-solitons meet at the origin at f = 0. 

For N - 3, any irreducible and totally non-negative A-matrix has the form. 



A = 



10-/7 

1 fl 



where a, b and c are positive constants. The r-function is then given by 

T = £(l,2,3) + fl£(l,2,4)-HZ7£(l,3,4)-Hc£(2,3,4), 

with E(l,m,n) - (A;„ - km){kn - ki)(km - ki)EiE„E„. This gives (1, 3)-soliton solution which is dual to the case of 
N - 1, that is, [l,4]-soliton fory <k and [1,2]-, [2,3]- and [3,4]-solitons fory » 0. The corresponding label for 
this solution is given by tt = (2 3 4 1). The lower figures in Figure 3.2 shows the time-evolution of this type. Here we 
set a, b, c of the A-matrix as 

|1,2,3| , |1,2,3| |1,2,3| 



a - 



b = 



11,2,41' 11,3,41' 12,3,41' 

where \l, m, n\ :- \(k„ - k„)(k„ - ki)(k^ - ki)\, so that all four solitons meet at the origin at f = 0. 



3.3. N — 2 and M — 4 cases 

By a direct construction of the derangements of ^4 with two exedances (i.e. = 2), one can easily see that there 
are seven cases with 2x4 irreducible and totally non-negative A-matrices. Then the classification theorems imply 
that we have a (2, 2)-soliton solution associated to each of those A-matrices. We here discuss all of those (2, 2)-soliton 
solutions with the same A; -parameters given by {k\,k2, kj,,kn) - 



■|, |, j), and show how each A-matrix chooses a 
particular set of line-solitons. In Figure [33] we illustrate the corresponding chord diagrams for all those seven cases. 



from which one can find the asymptotic line-solitons in each case. 



3.3.1. The case 7T — (3412) 

From the chord diagram in Figure [33] one can see that the asymptotic line-solitons are given by [1,3]- and [2,4]- 
types for both y — > ±00. The solution in this case corresponds to the top cell of Gr(2,4), and the A-matrix is given 
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(3412) 





(4321) 





(2413) 



(3142) 



(2143) 



Figure 3.3: The chord diagrams for seven different types of (2, 2)-soliton solutions. Each diagram parametrizes a unique cell in the totally non- 
negative Grassmannian Gr(2, 4). 



by 



A = 



I -c -d 
Q \ a b 



where a,b,c,d > are free parameters with D :- ad - be > 0. This is the generic solution on the maximal 
dimensional cell, and it is called T-type (after |9|). The most important feature of this solution is the generation of 
four intermediate soli tons forming a box at the intersection point. Those intermediate solitons are identified as [1,2]-, 
[2, 3]-, [1,4]- and [3, 4]-solitons, and they may be obtained by cutting the chord diagram of (3412)-type at the crossing 
points. For example, if we cut the chords of [1, 3] and [2,4] at the crossing point, we obtain either pair of {[1, 2], [3,4]) 
or {[1, 4], [2, 3]). The first pair appears in the lower and upper edge of the box, and the second pair in the right and left 



edges. Figure 3.4 shows the time-evolution of the solution of this type. The parameters a, b, c and d in the A-matrix 
determine the locations of solitons, their phase shifts and the on-set of the box-shaped interaction pattern; Those are 
given by 



|2,1| 0. 
14, 11 



|2,1| 

13, 2r 



0+ 



a 
d 



|4,2| 
I3,l|' 



D 



|2,1| 
|4,3| 



(3.3) 



where s is for the phase shift, j-^ for the location of [;, y]-soliton for x » (recall the definition below (2.10 1, i.e. 

Vj] ~ ~\h®l'J] S'^^^ x-intercept of the line of the [i, 7]-soliton for x » at f = 0), and r for the on-set of the 

1 



box. The phase shifts for [1, 3]- and [2,4]-solitons are given by Ax[,j] : = 

|1,4||2,3| be] 



VJ] 



VJ] 



\jJ\ 



©T with 



©T = In 



1,2||3,4| D 



(3.4) 



In Figure 3.4 we have chosen those parameters as s - 



1.3] 



forms an X-shape without phase shifts and opening of a box at the origin. In Figure 



® [2 4] ~ ^ ^'^'^ ^ ~ ^ ' ^^'^^ at f = the solution 
the amplitudes are A[i3] = ^ 



and A[2,4] - ||. Among the soliton solutions for N - 2 and M - 4, this solution is the most complicated and 



3.4 




40 -20 20 40 



-40 -20 20 40 



Figure 3.4: Example of (3412)-type soliton solution (T-type). The asymptotic solitons are [1,3]- and [2, 4]-types for both y 
intermediate solitons forming a box-shape pattern are given by [1, 2]-,[l, 4]-, [2, 3]- and [3,4]-solitons. 



±oo. The 



interesting one. As you can see below, this solution contains all other solutions as some parts of this solution, that is, 
as explained above, all six possible solitons (i.e. = 6) appear in this solution. 
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3.3.2. The case 7T ^ (4312) 

From the chord diagram in Figure 3.3 one can see that the asymptotic Hne-soUtons are given by [1,4]- and [2, 3]- 
solitons for y » and [1,3]- and [2, 4]-solitons in y <K 0. The A-matrix in this case is given by 



I -b -c 
1 fl 



where a,b,c > Q are free parameters. Figure 3.5 illustrates an example of this solution. The parameters in the A-matrix 




Figure 3.5: Example of (4312)-type soliton solution. The asymptotic solitons are [1,4]- and [2, 3]-types for y 2> and [1,3]- and [2, 4]-types for 
y <K 0. The intermediate solitons are [3, 4]-type for t = -6 and [1, 2]-type for t = 6. 

determine the locations of the line-solitons in the following form, 



12' 11 0.3, 
a — e f"i, 

13, 11 



^ |2,1| 0- ^ 13, 1| 0. 
b - — — e ['.31 and c = ---r-e I' -'i. 



13, 2| 



|4,3| 



The Q^.^j represent the locations of the [/, j] -soliton for x » and x «c 0, and the three line-solitons determine the 
location of the other one. We here take the parameters = 0|*j = - 0, so that all those four solitons meet 
at the origin at f = 0. Notice here that the pattern in the lower part (j < 0) of the solution is the same as that of the 
previous one, T-type. This can be also seen by comparing the chord diagrams of those two cases. 

5.5.5. The case 7T — (3421) 

The asymptotic line-solitons in this case are given by [1,3]- and [2,4]-solitons fory » and [1,4]- and [2,3]- 
solitons in 37 <«c 0. The A-matrix has the form. 



A = 



1 -c 
I a b 



where a,b,c > are free parameters. Figure 3.6 illustrates the evolution of the solution of this type. The parameters 




Figure 3.6: Example of (3421)-type soliton solution. The asymptotic solitons are given by [l,3]-and [2, 4]-types for y 2> and [1,4]- and 
[2, 3]-types for y <K 0. This solution is dual to the previous one of (4312)-type, i.e. the duality is given by (x,y, t) — > (-,v, -y, -f). 



in the A-matrix are related to the locations of the line-solitons. 



|2,1| ^ |2,1| e- ^ |2,1| 0- 
= — -e pji, b - ——e^M and c- e I'-'i. 



13, 1| 



|4,1| 



|4,2| 
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This solution can be considered as a dual of the previous case (b), that is, two sets of line-solitons for y » and y <K 
are exchanged. Here we take ©j^ - for those pairs (/, j) - (1 , 4), (2, 3) and (2, 4), hence all those four solitons meet 
at the origin at f = 0. Note here that now the upper pattern of the solution is the same as that of the T-type. So if one 



cuts Figures 3.5 and 3.6 along the x-axis and glues the lower half of Figure 3.5 together with the upper half of Figure 



3.6 we can obtain Figure 3.4 (T-type). This can be observed from the corresponding chord diagrams. 



One should also note that this solution is "dual" to the previous one in the sense that the patterns of those solutions 
are symmetric with respect to the change ix,y, t) (-x, -y, -t). 



3.3.4. The case n ^ (2413) 

The asymptotic line-solitons are given by [1,2]- and [2,4]-solitons for y » and [1,3]- and [3,4]-solitons for 
y <sc 0. The A-matrix is given by 

/I -c -d\ 
1^0 1 a bj' 

where a,b,c,d > with ad - be = 0. Those parameters are related to the locations of the Une-solitons, 



bJ^e%... 



0=7^^6*^11.31 and - = i^l^e®n,2i. 



d _ |4, 1| 
b ~ \4a\' 



|4,1| |3,2| 

with a = bc/d. Figure fiJ\ illustrates the evolution of the solution of this type. For the case that all the solitons meet at 




Figure 3.7: Example of (2413)-type soliton solution. The a.symptotic solitons are [1,2]- and [2, 4]-types for )' 2> and [1,3]- and [3,4]-types for 
)• <K 0. The intermediate solitons are [1, 4]-type in t = -6 and [2, 3]-type int = 6. 



the origin at f = 0, we set all to be zero. Notice that the nonlinear interaction generates the intermediate soliton 
of [1,4] -type in f < 0. This has the maximum amplitude among the (local) solitons appearing in the solution. It is 
interesting to note that the maximum amplitude soliton might be expected from the chord diagram, that is, this soliton 
is generated by the resonant interactions of [ 1 , 2] -H [2, 4] = [ 1 , 4] for y > and [1,3] + [3,4] = [1,4] for y < 0. One 



should also note that the left patterns in Figure 3.7 including [l,4]-soliton for f < and [2, 3]-soliton for f > are the 



same as those in Figure 3.4 (T-type). 



3.3.5. The case 7T — (3142) 

The asymptotic line-solitons are given by [1,3]- and [3, 4] -solitons in y » and [1,2]- and [2,4]-solitons in 
y <K 0. The A-matrix has the form, 

(1 a -c\ 
^"^0 1 bj' 

where a,b,c > 0. We represent those parameters in the following form which is particularly useful for our stability 
problem discussed in the next section. 



13, 1| 0. 

se r-A] 

|3,2| 



13, 1| 0^ 
.se 11-31 

14, 11 



and 



13, 1| 
|4,3| 



(3.5) 



Here represent the location of the [/, j] -soliton in x » 0, and s represents the location of the other Une-solitons 
in X <K 0, that is the generation of the intermediate line-solitons either [1,4]- or [2, 3]-type (c.f. the case (c) in Section 
|5]l. This solution is dual to the previous one of tt = (2413). In fact. Figure [378] illustrates the dual behavior of the 
evolution in the sense of the reverse of time. Here we take s = 1 and = ®12M ~ ^° ^^'^^ solitons meet 
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at the origin at f = 0. As in the previous case, the large amplitude soliton of [1,4] -type is generated by the resonant 
interaction. In particular, this solution can be used to explain the Mach reflection observed in shallow water ifTTl lSllSl. 




3.3.6. The case n — (4321) 

The asymptotic hne-solitons are [1,4]- and [2, 3]-types for both y — > ±oo. Since those solitons can be placed both 
in nearly parallel to the y-axis, this type of solutions fits better in the physical assumption for the derivation of the KP 
equation, i.e. a quasi-two dimensionality. This is then referred to as P-type (after |i9J). The situation is similar to the 
KdV case, for example, two solitons must have different amplitudes, A[i 4] > A[2,3]. The A-matrix is given by 



10 0-^ 
1 fl 



The parameters are expressed by 



a = ^ ^e'='M and b = ^e®i..4i. 



|4,3| 



|4,3| 



(3.6) 



where ©j^^.j represents the location of the [/, j] -soliton in x > (the wave front) (notice the phase shifts for the 
solitons). Figure '. 



3.9 



illustrates the evolution of the solution of this type. We here set - 0, and hence the solitons 




Figure 3.9: Example of (4321)-type soliton solution (P-type). The asymptotic solitons are [1, 4]- and [2, 3]-types for both y — > ±00. The pattern of 
the solution is stationary, and the interaction generates a negative phase shift of the solitons. 

in jc > (the wave front) meet at the origin at f = 0. Note here that there is a phase shift due to the interaction, and the 
shift is in the negative direction due to the repulsive force in the interaction similar to the case of KdV solitons. The 
phase shifts are given by Ax[, j] - -jjj^Qp with 



©P = In 



|4,2||3,1| 
|2,1||4,3| 



> 0. 



(3.7) 



This negative phase shift is due to the repulsive force in the interaction as in the case of the KdV solitons. One 
should also note that the interaction point moves in the negative y-direction, because the larger soliton of [l,4]-type 
propagates faster than the other one. Those [1,4]- and [2, 3] -solitons only appear as the intermediate solitons in the 



solution of T-type, that is, they form a part of the box in Figure 3.4 
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3.3.7. The case 7T ^ (2143) 

The asymptotic line-solitons are [1,2]- and [3,4]-types for both y ±oo. This solution was used originally 
to describe two soliton solution, and it is called O-type (after |9|). Interaction properties for soli tons with equal 
amplitude has been discussed using this solution. However this solution becomes singular, when those solitons are 
almost parallel to each other and close to the y-axis, contrary to the assumption of the quasi-two dimensionality for 
the KP equation |J IJ . The A-matrix is given by 



1 fl 

1/7 



Notice that this A-matrix is the limit of that for (3142)-type with c — » and k2 — » ^-3. We will further discuss this 
issue in the next section when we study the initial value problem with V-shape initial wave. The parameters a and b 
are expressed by 

|4,1| 0. 
a - --— e I' -i 



|4,2| 



and /,= Me0M. 



|4,2| 



(3.8) 



where ®^[i,j] represents the location of the wave front consisting of [1,2]- and [3,4]-solitons. Figure 3.10 illustrates 
the evolution of this type solution. Here we take all 0|^. j-^ = 0, the solitons in jc > (the wave front) meet at the origin 




Figure 3.10: Example of (2143)-type soliton solution (O-type). The asymptotic solitons are [1,2]- and [3, 4]-types for both y — » ±00. The pattern 
of the solution is stationary, and has a positive phase shift in the solitons. 

at f = 0. The phase shifts of those solitons are positive due to an attractive force in the interaction, and they are given 
by Axyj] = -1^100 with 

/|4, 1||3,2|\ 

When ^2 approaches to k^, the phase shift gets larger, and the interaction part of the solution is getting to be close to 
[l,4]-soliton as a result of the resonant interaction with [1,2]- and [3,4]-solitons at ^2 = ^3- This is the main argument 
of the paper till about the discussion on the Mach reflection of shallow water waves. Note also that those [1,2]- and 



[3, 4] -solitons appear as the intermediate solitons forming a box of T-type as shown in Figure 3.4 



4. Numerical simulations 

The main purpose of the numerical simulation is to study the interaction properties of line-solitons, and we will 
show that the solution of the initial value problem with certain types of initial waves asymptotically approaches to 
some of the exact solutions discussed in the previous section. This implies a stability of those exact solutions under 
the influence of certain deformations (notice that the deformations in our initial waves are not so small). 

The initial value problem considered here is essentially an infinite energy problem in the sense that each line- 
soliton in the initial wave is supported asymptotically in either y » or y <K 0, and the interactions occur only in 
a finite domain in the xy-plane. In the numerical scheme, we consider the rectangular domain D = {ix,y) : \x\ < 
L^, \y\ < Ly), and each line-soliton is matched with a KdV soliton at the boundaries y = ±Ly. In this section, we 
explain the details of our numerical scheme and give an error estimate of the scheme. 
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4.1. Numerical scheme 

Using the so-called window technique IfTTl [T9]| . we first transform our non-periodic solution u having non- 
vanishing boundary values into a function v rapidly decaying near the boundaries y - +L,, so that we consider the 
function v to be periodic in both x and y. With the window function Wiy), this can be expressed by the decomposition 
of the solution u, i.e. 

M = V + (1 - W)u with V - Wu. 



We take the window function Wiy) in the form of a super-Gaussian, 



W{y) = exp -Qn 



(4.1) 



where a„ and n are positive constants (we choose a„ - (1.1 1 1)" In 10 with n -21 fW\). 

We assume that near the boundaries y - +Ly, the solution can be described by exact line-solitons (i.e. we take Ly 
large enough so that the interaction will not influence those parts near boundaries up to certain finite time). We then 
consider the transformation in the form, 

u = v + {\-W)uQ, (4.2) 
where mq consists of exact line-solitons satisfying 

(1 - W)uo = (1 - W)u. 

This assumption prevents any disturbance from the boundary, and keeps the stability of the far fields consisting of 
well separated line-solitons. We are interested in studying an asymptotic behavior near the interaction point, and we 
do not expect to see any strong effect from the boundary region in a short time (recall that the system has a finite group 
velocity for waves with non-zero A:^ -component). 

With the transformation (|4.2|i, we have the equation for v. 



(45,v + 6vd^v + dlv) + 3d]v 
= 6(1 - W)d, (Wmo5,.mo - 5,.(vMo)) + 3(2W'dyUo + H^"mo). (4.3) 



Assuming the solution v to be periodic in both x and y with zero boundary condition at \y\ — L,,, we solve ( |4.3| l by 
using fast Fourier transform (FFT). For a convenience, let us rescale the domain D - {{x,y) : \x\ < L^, \y\ < Ly} with 
D := {(X, Y) : \X\ < n, \Y\ < n], i.e. 

n n 
X = —X, Y = —y. 

'-'X '-'y 

Then ( |4.3[ ) becomes 

dx (<9,v + Pvdxv + Qdlv) + Rd\v = F, (4.4) 

where F is the right hand side of ( |4.3| l with the scalings, and 

37: ^ TV' , 3;rLr 

P^Wl-^ G = 7TT' and 7? = 



^ 4Lr 4L2 ■ 

Under the periodic assumption of v, we have 



oo oo 

J(lX+mY) 



v{x, y, f) = ^ 



/=-oo m=-o 

Here the Fourier transformation v :- Tiv) is given by 



v{l, m, t) = f [ v(X, Y, t)e^''^'^^'"^^dXdY 

(2kY J_„ J_„ 
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Figure 4.1: Numerical accuracy test on T-type solution shown in Fig.|3.4| The graph shows the estimate of max|?^(v') - as a function of 

the number of Fourie modes N X N/2 for t = 0,3 and 6. The error is decreasing exponentially with respect to N until it saturates at round of error 
a 10"'* for N a= 1280, 1472 and 1536 for t = 0,3 and 6, respectively. The numerical domain is D = [-64, 64] X [-32, 32]. 



and the equation becomes 
with 

Let c = ^ - QP and V = e'"v. Then we have 

d,V + il-e'"N(v) = -e'"riF). (4.5) 
2 // 

We write this equation in the form of time evolution, 

d,V = Git, V). 

and use the 4-th order Runge-Kutta (RK4) method to solve the equation: Namely denoting f„ = nAf and 
Ki^G (f„, y„) , K2^G(t„ + ^Af, y„ + I^Af j , 
Ki^G (t„ + ^ Af, y„ + Af j , 7i:4 = G (tn + At, V„ + Y Ar) , 

we have 

= V« + z [^1 + 2(^2 + ^3) + ^^4] Af. 



We use a pseudo-spectral method to evaluate nonlinear term N(v) in ( |4.5| l and the 2nd order central difference to 
approximate three derivative terms in the F on the right hand side of ( |4.4| i for simplicity. This is because uq in ( |4.2| i 
is a non-periodic function and the estimation of derivatives from a Fourier spectral method still introduces errors 
at the computational boundary. A better choice is to provide the exact closed forms for the jc-and y-derivatives of 
Mo- Since some of the exact solutions have rather complicated formulae, we simply use 2nd order central difference 
approximation. Thus the numerical solutions are expected to be of 2nd order accuracy when there are enough Fourier 
modes in space and the solutions are vanished in large x. 

In order to verify the accuracy of the method, we provide two tests. The first test intends to answer the question on 
how many Fourier modes are needed to ensure machine accuracy in space for the solution v of (|4.3[). The second test 



shows the overall algorithm is 2nd order accurate. In Figure 4. 1 we show the maximum error in the Fourier modes of 



the T-type solution in the sense that we estimate max|;F(v) - ;F(v'exact)l with respect to the number of Fourier modes 



in both X- and y-directions for the case shown in Figure 3.4 On the computational domain D - [-64, 64] x [-32, 32], 
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Af = 10-^ 


Af = 5 X 10-^' 


order 


t = 1 


1.046 X 10-2 


2.645 X 10-^' 


1.98 


t = 2 


1.893 X 10-2 


4.353 X 10-^ 


2.12 


f = 3 


2.507 X 10-2 


5.986 X 10-^ 


2.07 


/ = 4 


2.992 X 10-2 


7.732 X 10-^' 


1.95 


t = 5 


3.523 X 10-2 


9.675 X 10-^' 


1.86 


f = 6 


4.211 X 10-2 


1.194 X 10-2 


1.82 



Table 1: The convergence test on the domain D = [-128, 128] X [-16, 16] with 2048 X 256 Fourier modes. The graph shows the error \\u' - 
"cxactllL2(D)/ll"exactllL-(0) ^™ Steps Af = 0.01 and 5 X IQ-'. This indicates that our scheme is of 2nd order accurate in time. 



machine accuracy « 10-'^ is reached for N ^ 1280, 1472 and 1536 for f = 0, 3 and 6, respectively. Since the T-type 
solution develops a box at later time, it requires more Fourier modes to well represent the solution. From this test 
with = 1536, we find that we need roughly = 12 Fourier modes per 1 unit in x to reach 10-'^ accuracy. Notice 
that it might take longer time to run such a fine mesh for a very large computational domain. If we choose to use 8 
Fourier modes per unit length in both x- and y-directions, the spacial error is within 10-'° and it can also speed up 
the computation. In Table [T] we summarize the error of the solution given by \\u' - Mexactllz.2(D)/l|MexactllL-(D) f^i" '^^^ 
different time steps Af = 0.01 and Af = 5 x 10-^ on D = [-128, 128] x [-16, 16] with 2048 x 256 Fourier modes. It 
can be clearly seen that the order of convergence is 2 because the error decreases by a factor around 4 when the time 
step is halved. If we take Af = 5 x 10-^ and keep running the simulation up to f = 12, the magnitude of the error is 
always within 2.5 x lO ^. In the following simulations, the errors of the solutions are at this magnitude or smaller. 



5. Numerical results 



Now we present the results of the direct numerical simulation of the KP equation with the following two types of 
initial waves consisting of two line-solitons. In particular, we fix the amplitude of one of the line-solitons, and place 
these line-sohtons in the symmetric shape with respect to the x-axis. This then implies that each profile of the initial 
wave is determined by two parameters Aq (the amplitude of other soliton) and (the angle of the wave-vector of the 
soliton): 

(i) V-shape wave consisting of two semi-infinite line-solitons, i.e. u(x, y, 0) = m j (x, y) + u^{x, y) with 

(An 



ii\{x, y) - Aq sech" 



-(x-ytanl'o) x H(y), 



(5.1) 



M°(x, y) = 2 sech2(x + y tan ^Pq) x H{-y), 



where Hiz) is the unit step function. 



if 
if 



z>0, 
z < 0. 



Here one should note that the initial locations x^^ ^.j are zero for both line-solitons, that is, they cross the origin 
at f = 0. The shifts of x^^. observed in the simulations provide the information of the interaction property. 

(ii) X-shape wave of the linear combination of two line-solitons, i.e. u{x,y, 0) = u^y(x,y) + u^{x,y) with 

[An 



M°(x,y) = Ao sech^ -J—(x-yian^o), 
ui{x, y) -2 sech2(x + y tan ^Fq). 



(5.2) 



Again note that the initial locations are x^. ., = 0. 
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These initial waves have been considered in certain physical models of interacting waves [[TTl |5] [161 HI El- In 
particular, in lfT6l[T8l . the V-shape initial waves were considered to study the generation of large amplitude waves in 
shallow water 

The main result of this section is to show numerically that those initial waves converge asymptotically to some of 
the exact solutions presented in Section [3] Here what we mean by the convergence is as follows: We first show that 
the solution u'(x, y) at the intersection point of V- or X-shape generates a pattern due to the resonant interaction, and 
we then identify the pattern near the resonant point with that given by some of the exact solutions u'^^^^^(x,y). The 
convergence is only in a local sense, and we use the following (relative) error estimate. 



E{t) \\U' - «Lactlli:(fl^)/ll«exactll22(B,), (5.3) 

where DJ. is the circular disc defined by 

d;. := {(x,y) eR^:(x- xo(t)f + (y - yo(t)f < r^) . 

Here the center (xo(f),3'o(0) is chosen as the intersection point of two lines determined from the corresponding exact 
solution identified by the following steps: 

1. For a given initial wave, we first identify the type of exact solution (i.e. the chord diagram) from the formulae, 
A[, y] = j(ki - kj)^ and tan^^ j] = k, + kj. This then gives the exact solution with the corresponding A-matrix. 
Note here that the entries of the A-matrix are not determined yet. 



2. We then minimize the error function E{t) of ( |5.3| l at certain large time f = by changing the entries in the 
A-matrix, that is, we consider 

Minimize £(70). (5.4) 

A-matrix 

The minimization is achieved by adjusting the solution pattern with the exact one, that is, the pattern determines 
the A-matrix whose entries give the locations of all line-solitons including the intermediate solitons in the 
solution. With this A-matrix, the center (xoit),yoit)) of the circular domain D',. is given by the intersection point 
of the lines, 

x + ytan *l'[i,j,] - C [,-, j,]f = jc+ for / = 1 , 2, (5.5) 

where the locations x'^^j^ j^-^ are expressed in terms of the A-matrix (see Section [3|l. Here we have Aj,, j,] - 
A(), A[,, J,] = 2 and *l'[,,j,] - -^^[i, j,] = ^o- The time Tq may be taken as To - 0, out for most of the cases, we 
need To to be sufficiently large. This implies that the solutions have nontrivial location shifts, i.e. in 

the solution. 

3. We then confirm that E{t) further decreases for a larger time t > Tq up-to a time t - T\ > Tq, just before the 
effects of the boundary enter the disc DJ. (those effects include the periodic condition in x and a mismatch on 
the boundary patching). 

We take the radius r in D'. large enough so that the main interaction area is covered for all f < Tq, but DJ. should be 
kept away from the boundary to avoid any influence coming from the boundaries. The time > Tq gives an enough 
time to develop a pattern close to the corresponding exact solution, but it is also limited to avoid any disturbance from 
the boundaries for f < Ti . In this paper, we give To and Ti based on the observation of the numerical results, and we 
will report the detailed analysis in a future communication. 

5.7. V-shape initial waves 

The initial wave form is illustrated in the left panel of Figure 5.1 The right figure shows the (incomplete) chord 
diagrams corresponding to the initial V-shape waves; that is, the upper chord represents the (semi-infinite) line-soliton 
in y > and the lower one represents the (semi -infinite) line-soliton in y < 0. We have done the numerical simulations 



with six different cases which are marked in Figure 5. 1 as (a) through (f). The main result in this section is to show that 
each incomplete chord diagram representing the initial wave gets a completion with minimum additional chords which 
represents the asymptotic pattern of the solution. Here the minimal completion implies that for a given incomplete 
chord diagram, one adds extra chords with minimum total length so that the complete chord diagram represents a 
derangement (see also |10|). For our examples of the cases (a) through (f), we have the following completions: 
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Figure 5.1: V-shape initial waves. Each line of the V-shape is a semi-infinite line-soliton. We set those line-solitons to meet at the origin, and fix 
the amplitude of the soliton in y < to be 2. Aq is the amplitude of the soliton in y > 0. The right panel shows the incomplete chord diagrams 
corresponding to the choice of the values Aq and . There is no corresponding exact soliton solution on the dotted line of k2 = k^. 



(a) With k\ = ki, the completion is ;7r = (^^^), i.e. add [3,4]-chord iny > 0. This implies that [3,4]-soliton appears 
in the asymptotic solution in y>Q, and the corresponding exact solution is a (1, 2)-soliton (cf. Figure [TT) . 

(b) With k^ - k4, the completion is tt = (3 J 2), i.e. add [l,2]-chord in y < which appears in the asymptotic 



solution, and the corresponding solution is a (2, l)-soliton (cf. Figure 3.1 



(c) The completion is ;7r = (3142), i.e. add [l,2]-chord in y < and [3,4]-chord in y > 0. These corresponding 



solitons appear behind those initial solitons (cf. Figure 3.8 1 



(d) The completion is tt = (2143), i.e. add [1, 2]-chord in y < and [3, 4]-chord in y > 0. The soliton solution is of 



O-type, and the semi-infinite initial solitons extend to generate other half with the phase shift (cf. Figure 3.10i 



(e) The minimal completion is tt = (2341), i.e. add [1,2]- and [3,4]-chords in y > 0. This is a (1, 3)-soliton. (Note 
that the initial wave is a front half of P-type soliton solution, but the minimal completion is not of P-type.) Since 
[1, 2]-soliton iny > appears ahead of the wavefront consisting of [2, 3]-soliton iny > and [l,4]-soliton in 
y < 0, the generation of [1, 2]-soliton is impossible. The [1, 2]-soliton is a virtual one, but it gives an influence 
on bending the [l,4]-soliton in y > (see below for the details). 

(f) The minimal completion is ;r = (4123), i.e. add [1,2]- and [3,4]-chords in y < 0. The is a (3, l)-type, but again 
we note that [3, 4]-soliton in y < is a virtual one (i.e. it locates ahead of the wavefront similar to the previous 
case (e)). We also observe the bending of [l,4]-soliton in y > under the influence of this virtual soliton (see 
again below for the details). 



5.7.7. The case (a) 

This is a critical case with ki = k2. The line-solitons of V-shape initial wave are [1, 3]-soliton in y > and 
[ 1 , 4] -soliton in y < 0. We take the following values for the amplitudes and angles of the those solitons, 

A[i.3] = Ao = 5, tan ^[1,3] = = -tan>Po, 
A[i,4]=2, tan»P[i,4] = i = tan^o. 

With the formulas for the amplitude A[,j] - ^(k, - kj)^ and the angle tan*F[, j.] - k, + kj, we have the A:-parameters 
{k\ - k2, A;3, A;4) = (-|, |, |). The corresponding exact solution is of (1, 2)-type. 



The top figures of Figure 5.2 illustrate the result of the direct numerical simulation. The simulation shows a 
generation of third wave at the point of the intersection of those initial solitons. This generation is due to the resonant 
interaction, and this may be identified as [3,4]-soliton of the exact solution, that is, the resonance condition gives 
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Figure 5.2: Numerical simulation of the case (a) in Figure [TT] The initial wave consists of [1, 3]-soliton in y > and [1, 4]-soliton in y < 0. The 
top figures show the result of the dir ect si mulation, the middle figures show the corresponding exact solution of (1, 2)-type, and the bottom one 
represents the en'or function E{t) of \5.3\ which is minimized at f = 10. The circles in the exact solution shows the domain D[.. A dispersive 
radiation appearing behind the interaction region may be considered as a reaction of the generation of the third soliton of [3, 4]-type as a result of 
resonant interaction (see the text for the details). 



[1,3] + [3,4] = [1,4] and completes the incomplete chord diagram. The generation of [3,4]-soliton may be also 
understood by decomposing the exact solution into two parts, say Mexact - u + v, where u is the solution of the KP 
equation with our initial wave of [1, 3]- and [l,4]-solitons. Then v satisfies, 

5,-(45,v + 6vd,,v + dl-v) + 3d]v = -6dl(uv), 

which is the KP equation with one extra term in the right hand side. If we assume that the initial solitons remain as 
the part of the exact solution Mexact. the product uv in the extra term has nonzero value only near the origin at f = 
(this assumption is not totally correct, but the shifts of the initial solitons are very small for most of the cases). We 
then expect v to be close to the solution of the KP equation with the initial wave of [3, 4] -soliton only for y > Q and 
X < 0. Since this [3,4]-soliton is missing in the front half {x > 0), the edge of the wave gets a strong dispersive effect. 
Then the edge of the solution v disperses away with a bow-shape wake propagating toward the negative x-axis as can 



be seen in the simulation (cf Figure 5.2 1. Recall here that the dispersive waves (i.e. non-soliton parts) propagate in 
the negative x-direction (see Section |3]l. Then the decay of the v solution at the front edge implies the appearance of 
the [3,4]-soliton in the exact solution as one can see from u - Mgxact - v. This argument may be also applied to other 
cases as well. 

Minimizing the error function E{t) at f = To - 10 with ( |5.4| i, we obtain the A-matrix of (1, 2)-type, 

/I -1.27 
1^0 1 0.63 

The minimization is achieved by adjusting the locations of the [1, 3]- and [1,4] -soli tons of the wavefront, and it gives 
the shifts 

x+ 3j = -0.0079, x+ 4j = -0.120. 
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Figures. 3: Numerical simulation of the case (b) in Figure |5.1| The initial wave consists of [1, 4]- and [2, 4]-solitons in v > Oandy < respectively. 
The top figures show the result of the direct simulation, and the middle figures show the corresponding exact solution of (2, l)-type. The bottom 
graph shows the en'or function E{t) of j53j which is minimized at ; = 6. 



(The A-matrix is then obtained using (3.2i). Thus the new Hne-soHton of [3,4]-type has a relatively large negative 
shift. This may be explained as follows: The [3, 4]-soliton is generated as a part of the tail of [1, 4]-soliton, and in the 
beginning stage of the generation, the [3, 4] -soli ton has a small amplitude and propagates with slow speed. This causes 
a negative shift of this soliton. Also the negative shift in the [l,4]-soliton is due to the generation of the [3, 4] -soli ton, 
i.e. losing its momentum. The middle figures of Figure 5.2 show the exact solutions generated by the r-function with 



those parameters (^i = k^, k^, ki^) - (-|, \, |) and the A-matrix given above. The circles in the figures show D\ where 
we take r = 12 and the center is given by (-Xo(0 = ~ 0.0637, yoCO - \^ ~ 0.1 12). In the bottom figures of Figure 
5.2 we plot the error function E(t) of (|5.3|l. Figure 5.2 indicates that the asymptotic solution seems to converge to the 



exact solution of (1, 2)-type with the above A-matrix and the same ^-parameters. 
5.1.2. The case (b) 

This is also a critical case with kj - k^,. The line-solitons of V-shape initial wave are [1, 3]-soliton in y > and 
[2, 3]-soliton in y < 0. We take the amplitudes and angles of those solitons to be 

I A[i,3] = Ao = f , tan'Pii^j = -| = -tan*Po, 
1 A[2,3] =2, tan>I'[2,3] = | ^tan*!',,. 

illustrates the result of the numerical 



5.3 



This then gives the A:-parameters (^1,^2, ^3 - ^4) - (~T'~8' T^- Figure 
simulation. The top figures show the direct simulation of the KP equation. Notice again that a bow-shape wake behind 
the interaction point corresponds to the v solution in the decomposition, u - Mexact - v, and the decay of v implies 
the appearance (or generation) of [1, 2]-soliton. Notice that v gives the [1, 2]-soliton part in the exact solution. The 
middle figures show the corresponding exact solution of (2, l)-type, whose A-matrix is obtained by minimizing the 
error function E{t) at f = 6, 

A = (l 1.15 1.21). 
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Figure 5.4: Numerical simulation of the case (c) in Figure \5A\ The initial wave consists of [1, 3]-soliton in y > and [2,4]-soliton in y < 0. 
The top figures show the result of the direct simulation, and the middle figures show the corresponding exact solution of (3142)-type. The bottom 
graph of the en'or function E(t) which is minimized at f = 10. A large amplitude intermediate soliton is generated at the intersection point, and it 
corresponds to the [ 1 , 4]-soliton with the amplitude A[i_4] = 4.5. 



The shifts of locations x^. obtained in the minimization are given by 



[1,3] 



;c+ 3j = -0.0254, 



using (3.1 1. Those negative shifts may be explained in the similar way as in the previous case. In those figures, the 
domain D[ has r = 12 and its center (xo(f) = - 0.0399, yo(t) = -yf + 0.0194). The bottom graph in Figure 
shows that the solution converges asymptotically to the exact solution of (2, l)-type with the above A-matrix and 
same A:-parameters. 



5.3 
tEe 



5.1.3. The case (c) 

The line-solitons of the initial wave are [1, 3]-soliton in y > and [2,4]-soliton in y < 0. We consider the case 
with the larger soliton in y > 0, and take the amphtudes and angles of those solitons as 

f A[i,3] = Ao =2, tan ^[1,3] = -1 = -tan'Po, 
\ A[2,4] = 2, tan ^[2,4] = 1 = tan>Po- 

The A;-parameters are then given by (ki,k2,k2,k4) - (-1, -I, I, I). The exact solution obtained by the completion of 



the chord diagram is given by (3142)-type soliton solution. Figure 5.4 illustrates the result of the numerical simulation. 
The top figures show the direct simulation of the KP equation. We observe a bow-shape wake behind the interaction 
point. The wake expands and decays, and then we see the appearance of new solitons which form resonant interactions 
with the initial solitons. One should note that the solution generates a large amplitude intermediate soliton at the 
interaction point, and this soliton is identified as [l,4]-soliton with the amplitude A[i,4] = 4.5. This solution has been 
considered to describe the Mach reflection of shallow water waves, and the [1,4] -soliton describes the wave called the 
Mach stem ||3l[T0l. 
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The middle figures in Figure 5.4 show the corresponding exact solution of (3142)-type whose A-matrix is found 
by minimizing E{t) at f = 10, 

h 1.92 -L96\ 
1 0.64 



A = 



In the minimization to find the A-matrix, we adjust the locations of solitons including the intermediate solitons, and 
using p.5|l, we obtain 



a;+ 3j = jt:+4j = -0.01, s = 0.980. 

Those values indicate that the solution is very close to the exact solution for all the time. The negative value of 
the shifts ^.j is due to the generation of a large amplitude soliton [l,4]-type (i.e. the initial solitons loose their 
momentum), and .? < 1 implies that the [1,4] -soli ton is generated after f = 0. Also note that the [1,4] -soli ton now 
resonantly interact with [1, 3]- and [2,4]-solitons to create new solitons [1,2]- and [3, 4] -solitons (called the reflected 
waves as considered in the Mach reflection problem lflTl [3l [T0l ). This process then seems to compensate the shifts 
of incident waves, even though we observe a large wake behind the interaction point. Notice that the wake disperses, 
and the exact solution starts to appear as the wake decays (recall again the decomposition Mexact = u + v (where v is 
the wake and u is the solution). 

The bottom graph in Figure [54| shows a rapid convergence of the initial wave to (3142)-type soliton solution with 
those parameters of the A-matrix and k values given above, and for the error function, we use D'^ with r = 12. 

5.7.4. The case (d) 

The line-solitons of the initial wave are [1, 2]-soliton in y > and [3,4]-soliton iny < 0. This case corresponds 
to the solitons with a large angle interaction. For simplicity, we consider a symmetric initial data with the initial wave 
whose amplitudes and angles are given by 

I A[i,2] = Ao = 2, tan»P[i,2] = = -tan^P,,, 
\ A[3,4] = 2, tan»l'[3,4] = y = tan^Po- 

The ^-parameters are then given by (ki,k2,k3,k4) = (~T'~I' 5' T-*- completion of the chord diagram gives 
O-type soliton solution, i.e. (2143)-type. This implies that the solution extends the initial solitons to the negative 
x-direction with the same types, i.e. [1,2]- and [3,4]-types in y < and y > 0, respectively. One should note that 



the solution also generates the phase shift determined uniquely by the ^-parameters. Figure 5.5 illustrates the result 
of the numerical simulation. The top figures show the direct simulation of the KP equation. The wake behind the 
interaction point has a large negative amplitude, and this corresponds to the v solution in Mexact = u + v. The middle 
figures show the corresponding O-type exact solution whose A-matrix is determined by minimizing the error function 
E(t) at f = 10, 

/I 1.91 \ 
1 0.17 



Because of our symmetric initial data, we consider the shifts x^^ - ^[j,4] for 'he minimization, and using (3.8 i, we 
obtain 

■"■[1,2] ~ -^[3,4] ~ -0.020. 

The negative shifts imply the slow-down of the incidence waves due to the generation of the solitons extending the 
initial solitons in the negative x-direction. Those new parts of soUtons have the locations x^i jy The phase shifts 
Ax[, j] = Xj^ ^.j - X|^ for the line-solitons are calculated from ( 3.9 1, and they are 



Ax[i,2] = Ax[3,4] = 0.593. 

The positivity of the phase shifts is due to the attractive force between the line-solitons, and this explains the slow- 



down of the initial solitons, i.e. the small negative shifts of xt .,. The bottom graph in Figure 5.5 shows E{t) which 



is minimized at f = 6, and we take r = 12 for the domain D'. One can see a rapid convergence of the solution to the 
O-type exact solution with those parameters. One should however remark that for a case near the critical one with 
k2 - kj,, there exists a large phase shift in the soliton solution, and the convergence is very slow. Note that the case 
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with k2 = ^3 also corresponds to the boundary of the previous case (c) (See Figure 5.1 1, and the amplitude of the 
intermediate soliton generated at the intersection point reaches four times larger than the initial solitons at this limit. 
This large amplitude wave generation has been considered as the Mach reflection problem of shaUow water wave 

5.1.5. The case (e) 

The line-solitons of the initial wave are [2, 3]-soliton in y > and [1, 4]-soliton in y < 0. We consider the solitons 
parallel to the y-axis with the amplitudes, 

I Ap,3] = Ao = 0.25, tan'Pp,3] = = - tan*Po, 
\ A[i,4] = 2, tan'P[i,4] = = tan*Fo- 

Then the ^-parameters are given by (^1,^2,^3, ^4) - ^1^° considered the case with nonzero 

angle (i.e. oblique to the y-axis with *Po > 0), and we expected to see a P-type soUton as the asymptotic solution. 
However, the numerical simulation suggests that we do not see a solution of P-type, but the asymptotic solution is 
somewhat close to the exact solution of (2341)-type as obtained by the minimal completion of the chord diagram. 
The (2341)-type is a (l,3)-type soliton solution with [l,4]-soliton in y < and [1,2]-, [2,3]- and [3, 4] -soli tons in 
y > 0. Here one should note that the [1, 2] -soliton in y > locates ahead of the wavefront consisting of [1,4]- and 
[2, 3]-solitons, so that this soliton cannot be generated by the resonant interaction of those [1,4]- and [2,4]-solitons. 



Figure 5.6 illustrates the result of the numerical simulation. The top figures show the direct simulation of the KP 
equation. We observe that the [1,4] -soliton in y < propagates faster and its tail extends to y > 0. Also note that 
a large dispersive wake appearing in y < has negative amplitude and disperses in the negative jc-direction. On the 
other hand, the tail part in y > seems to have a resonant interaction with the [2, 3] -initial soliton and becomes a new 
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Figure 5.6: Numerical simulation of the case (e) in Figure [TT] The initial wave consists of [2, 3]-soliton in y > and [1, 4]-soliton in y < 0. The 
upper figures show the result of the direct simulation, and the lower figures show the corresponding exact solution of (2341)-type obtained by the 
minimal completion of the chord diagram given by the initial wave. Notice that the bending on the initial [1, 4]-soliton appears near the point where 
the virtual [1, 2]-soliton intersects. The dispersive radiation behind the intersection point may be the part of this [1, 4]-soliton. 



soliton. This behavior of the solution may be explained by the chord diagram of (2341)-type whose corresponding 
exact solution is illustrated in the bottom figures in Figure 5.6 The new soliton generated in y > may be identified 
as [3, 4] -soliton. Also we note that the bending of the [l,4]-soliton observed in the simulation seems to appear at the 
point where the [1, 2] -soli ton has a resonant interaction with [1,4] -soli ton. This means that the intermediate soliton 
in the exact solution is [2, 4]-soliton, and this soliton may give a good approximation of the bending part of the [1,4]- 
soliton. It is interesting to consider this [1, 2]-soliton as a virtual one, which is not visible but has an influence to the 
initial solitons. This can be also understood by the similar argument using the decomposition Mexact - u + v. In this 
case, the initial data for v is given by the [1,2]- and [3,4]-solitons in Mexact- Then the [3,4]-soliton in Mexact appears 
when the solution v decays in this part. On the other hand, the decay of v in the [1, 2]-soliton part seems to occur only 
in the region behind the front solitons as seen at f = 5. Non-decaying of the [l,2]-soliton implies that [1, 2]-soliton 
never appears in the solution m, because of the cancellation of Mexact and v on this sohton part. We will give more 
details in a future communication. 



5.7.6. The case (f) 

The line-solitons of the initial wave are [1, 4]-soliton in y > and [2, 3]-soliton in y < 0. Contrary to the previous 
case, we have now the larger soliton in y > and the smaller one in y < 0. In a comparison, we also consider the 
solitons parallel to the y-axis. We then take the amplitudes and angles of the those sohtons as 

^[1,4] =^0=6, tan'I'fi^] = = -tan'I'o, 
A[2,3] = 2, tan 'I' [2,3] = = tan^Po. 

Then the A:-parameters are given Since the upper [1,4] -soliton propagates faster 

than the lower one, the lower tail of this soliton appears and then resonates with the lower [2, 3]-soliton. Following 
the similar arguments as in the previous case, we expect that (3, l)-soliton is the corresponding exact solution which 
is of (4123)-type as given by the minimal completion of the chord diagram. Figure [5?7| illustrates the result of the 
numerical simulation. The upper figures show the direct simulation, and the lower figures show the exact solution 
of (4123)-type. As in the previous case (e), this case shows that the [3,4]-soliton part in the exact solution does not 
appear in the solution of the simulation, that is, the [3,4]-soliton is a virtual one. This can be also explained with the 
similar argument as in the previous case. Notice again that the [1,4] -soli ton in y > has a virtual interaction with 
[3,4]-soliton to get bending near (x - l^,y - 7) at t - 6. 
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Figure 5.7: Numerical simulation of the case (f) in Figure [5?T] The initial wave consists of [1, 4]-soliton in y > and [2, 3]-soliton in y < 0. 
The upper figures show the result of the direct simulation, and the lower figures show the corresponding exact solution of (4123)-type obtained 
by the minimal completion of the chord diagram. Notice that the bending in the [1, 4]-soliton appears near the intersection point with the virtual 
[3, 4]-soliton, and the dispersive radiation observed in the simulation is a part of this [3, 4]-soliton. 




Figure 5.8: X-shape initial waves. Each line of the X-shape is an infinite line-soliton solution. We set those line-solitons to meet at the origin, 
and fix the amplitude of the soliton with positive angle To to be 2 (the corresponding chord is shown in the lighter color in each region). Aq is the 
amplitude of the soliton with negative angle -'Fq. The right panel shows the corresponding chord diagrams with the choice of the values Aq and 
Tq. There is no exact solution for the degenerate cases marked by the dotted lines. 



5.2. X-shape initial waves 

The initial wave form is illustrated in the left panel of Figure |5.8| The right figure shows the chord diagrams 
corresponding to the initial waves with the parameters (Ao,tan*Po)- We carried out four diff'erent cases marked with 



(g) through (j) in Figure 5.8 The main result in this section is to show that each simple sum of line-solitons gets 
asymptotically to the exact solution corresponding to the chord diagram given by the initial wave. One should however 
note that those initial waves are not the exact solutions, and the simulations provide stability analysis for those exact 
solutions. For example, if the exact solution has the phase shift, the asymptotic solution actually converges to this 
exact solution by generating the correct phase shifts. 



5.2.7. The case (g) 

The initial wave consists of [1, 2]- and [3, 4] -soli ton. To emphasize the phase shift appearing in the exact solution 
of (2143)-type, we consider the initial solitons to be close to the boundary at k2 - kj, (but we cannot take the values 
so close, since the generation of a large phase-shift needs a long time computation). For simplicity, we also take a 
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symmetric initial wave whose amplitudes and angles of the initial solitons are given by 



: Ao = 2, tan *i'[i,2j = -2.02 = - tan 
■ 1 tan »i'[3,4j =2.02 = tan >Po 

The A:-parameters are then given by {k\,k2,k-i„k\) - (-2.01,-0.01,0.01,2.01). The numerical simulation in Figure 




5.9 demonstrates the generation of the phase shifts after the interaction, that is, the line-soli tons in the left side are 
accelerated by their interaction. Note here that the interaction generates the dispersive wake behind the interaction 
point which propagates in the negative x-direction. 

Minimizing the error function E{t) at f = 12, we obtain the A-matrix, 



A = 



1 3.029 
1 0.015 



In the minimization by adjusting the locations of the solitons, we use x^^ - -^p 4] (due to the symmetry of the 
solution), and then we obtain 

01 — ^\-x A^ — —0.21. 



The negative values are due to the generation of the positive phase shifts given by ( 3.9 1, 

Ax[i_2] - Aji:[3,4] = 1.963 > 0. 

Namely the accelerations of the back parts of the line-solitons (i.e. the positive phase shifts) imply the deceleration 



of the front parts of solitons. The middle figures in Figure 5.9 shows the corresponding O-type soliton solution with 



this A-matrix. The circle shows the domain D\. with r - 12. Notice the dispersive wake behind the interaction point 



in the simulation. The bottom graph in Figure 5.9 shows the convergence of the solution to the O-type solution given 
by those A-matrix and the same A:-parameters. 

5.2.2. The case (h) 

The initial wave is the sum of [1, 3]- and [2, 4] -soliton. The corresponding chord diagram is of (3412)-type (i.e. 
T-type), and we expect the generation of a box at the intersection point of those solitons. For simplicity, we consider 
a symmetric initial wave, and take the amplitudes and angles of those solitons to be 



[1,3] 



= Ao=2, tan*P[i3] = -1 = -tan*Po 



A[2,4] = 2 tan*P[2,4] = 1 = tan*!',, 

The ^-parameters are then given by {k\,k2, k^,k^) - Although T-type soliton appears for smaller angle 

^Oj one should not take so small value. For an example of the symmetric case, if we take *Po - 0, we obtain the 
KdV 2-soliton solution with different amplitudes. So for the case with very small angle ^Pq, we expect to see those 
KdV solitons near the intersection point. However, the solitons expected from the chord diagram have almost the 
same amplitude as the incidence solitons for the case with a small angle. The detailed study also shows that near the 
intersection point for T-type solution at the time when all four solitons meet at this point (i.e X-shape), the solution 
has a small amplitude due to the repulsive force similar to the KdV solitons (see also |4|, where the initial X-shape 
wave with small angle generates a large soliton at the intersection point). This then implies that our initial wave given 
by the sum of two line-solitons creates a large dispersive perturbation at the intersection point (this can be also seen 
in the next cases where we discuss the P-type solutions). 



The top figures in Figure 5.10 illustrate the numerical simulation, which clearly shows an opening of a resonant 
box as expected by the chord diagram of T-type. The corresponding exact solution is illustrated in the middle figures, 
where the A-matrix of the solution is obtained by minimizing the error function E(t) of \5.3>\ at f = 6, 

ll -0.368 -0.330' 
io 1 1.198 0.123 
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Figure 5.9: Numerical simulation of the case (g) in Figure The initial wave is the sum of [1,2]- and [3,4]-solitons. The top figures show the 
numerical simulation, and the middle figures show the corresponding exact solution of (2143)-type, i.e. O-type. The bottom graph shows E(t) of 
{5.3\ which is minimized at ? = 12. Notice that a small dispersive radiation appears to generate the correct phase shifts. 



In the minimization, we take x; 



1,3] 



12,4] 



due to the symmetric profile of the solution, and adjust the on-set of the 



box (recall ( 3.3 i, and note that the symmetry reduces the number of free parameters to three). We obtain 



ni,3] 



'[2,4] 



0.025, 



3.63, 



s = 0.350. 



The positive shifts of those [1,3]- and [2,4]-solitons in the wavefront indicate also the positive shift of the newly 
generated soUton of [1,4] -type at the front. This is due to the repulsive force exists in the KdV type interaction 
explained above, that is, the interaction part in the initial wave has a larger amplitude than that of the exact solution, 
so that this part of the solution moves faster than that in the exact solution. This difference may result as a shift of the 



location of the [1,4] -soli ton. The phase shifts of the initial waves are obtained from (3.4i 



1,3] 



= Ajcp,4] = 0.549. 



The relatively large value r > 1 indicates that the on-set of the box is actually much earlier that f = 0, and s < I shows 
the positive phase shifts as calculated above. 

The bottom graph in Figure 5.10 shows the evolution of the error function E(t) of ( |5.3| l which is minimized at 
t - 6. Note here that the circular domain Dj. with r = 22 covers well the main feature of the interaction patterns for 
all the time computed for f < 7. 



5.2.3. The case (i) 

The initial wave consists of [1,4]- and [2, 3]-soliton. The corresponding diagram is of (4321)-type (P-type), and 
the exact solution of this type has negative (repulsive) phase shift. Here we show how the phase shift is generated in 
the simulation. We take the amplitudes and angles of those solitons to be 

A [2,3] = = ^, tan >P[2,3] = -0.4 = - tan ^Fq 
A[i_4] = 2 tan*F[i_4] = 0.4 - tan*Fo 
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Figure 5.10: Numerical simulation of the case (h) in Figure [5^ The initial wave is the sum of [1,3]- and [2, 4]-solitons. The top figures show the 
numerical simulation, and the middle figures show the coiTesponding exact solution of (3412)-type, i.e. T-type. The bottom graph shows the eiTor 
function E(t) of j53J which is minimized at ; = 6. 



The A:-parameters are then given by (fei, ^2,^3, ^4) = ~|' I' f )■ '^^is case, we also note that taking a very small 
angle ^'q may give a large perturbation to the exact solution, and one should have a large computation time to see the 



convergence. The top figures in Figure 5.11 show the numerical simulation. The simulation indicates the generation of 



the phase shift in the lower part of the solution. This may be explained as follows: Since the [1, 4]-soliton propagates 
faster than the [2, 3]-soliton, the intersection point moves in the upper direction, i.e. the positive y-direction. Then 
the effect of the nonlinear interaction may be observed in the lower half (i.e. behind the interactions) of the solution, 
and it generates the phase shifts. One should also note that the large bending of the [l,4]-soliton in y < 0. This is 
due to the interaction with the dispersive wake generated around the origin at ? = 0, and notice that this bending part 



disperses away. The middle figures in Figure 5.11 show the corresponding exact solution of P-type whose A-matrix is 
obtained by minimizing the error function E(t) of (|5.3[) at f = 12, 



A = 



1 
1 1.629 



-0.038 




In the minimization, we adjust the location of the line-solitons at the wavefront, and using p.6|l, we obtain 



^[2,3] 



= -0.125, 



4j — 0.54. 



Notice that the larger shifts in the [l,4]-soliton, which can be also observed in the simulation. This shift is due to the 
repulsive force appearing at the intersection point (recall that the exact solution has a smaller amplitude at this point). 



The phase shifts of the P-type solution is given by ( 3.7 1, 

Ax[i,4] = -1.10, Ax[2,3] = -2.75. 
Those large phase shifts are also observed in the simulation, and the bending regions are getting separated from the 



intersection point with the correct phase shifts. The bottom graph in Figure 5.11 shows the decay of the error function 
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Figure 5.11: Numerical simulation of the case (i) in Figure [5^ The initial wave is the sum of [1,4]- and [2, 3]-solitons. The top figures show the 
numerical simulation, and the middle figures show the corresponding exact solution of (4321)-type, i.e. P-type. The bottom graph shows the eiTor 
function E(t) of j53J which is minimized att = 12. Note that the bending on the solitons disperses and the interaction pattern converges locally to 
that of the corresponding exact solution. 



E(t) with D[ for r = 10. We also note the saturation E(t) after f = 1 1, which may be due to the large deformation of 
the solution, and we may need to perform a larger domain computation. 



5.2.4. The case (j) 

This case is similar to the previous one with [1,4]- and [2, 3]-soliton in the initial wave. However in this case, 
the larger soliton of [1,4] -type has a negative angle (propagating in the negative y-direction), and we expect that the 
phase shift now appear in the upper half of the solution. We take the amplitudes and angles of those solitons to be 




169 

32 ' 



tan*?. 



1.4] 



tan >P[2,3] 



= -tan*Po 
tan^Po 



5.12 



show the numerical 



The A:-parameters are then given by (ki,k2, k^, k4) = (-|, |, |). The top figures in Figure : 
simulation. One can observe a dispersive wake is generated by a large perturbation in the initial wave near the origin, 
and then the wake pushes the [1, 4]-soliton forward. Also the large bending of the [2, 3]-soliton near the origin is due to 
the interaction with the wake. Since the [1, 4]-soliton propagates much faster than the other one, the intersection point 
moves in the negative y-direction. We can see a clear separation between the dispersive waves causing the deformation 



of solitons and the intersection points of P-type. The middle figures in Figure 5.12 illustrate the corresponding P-type 
solution whose A-matrix is obtained by minimizing the error function E{t) at f = 8. We take the radius r - 10 for a 
good separation of the dispersive perturbation from the steady P-type soliton part. With this minimization, we get 



A = 



1 
1 8.549 



-0.299 




In the minimization, we adjust the locations of the line-solitons at the wavefront, and using (3.6 1, we obtain 



;c+_4j = 0.627, 



^[2,3] = 0.923. 
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Figure 5.12: Numerical simulation of the case (j) in Figure [5^ The initial wave is the sum of [1,4]- and [2, 3]-solitons. The top figures show the 
numerical simulation, and the middle figures show the corresponding exact solution of (4321)-type, i.e. P-type. The bottom graph illustrates the 
error function E(t) of \53\ which is minimized at f = 8. The graph shows a further reduction of the error, that is, the convergence of the solution 
to the corresponding exact solution even with a drastic change of the error at the initial stage. 



The positive shifts are due to the repulsive interaction at the intersection point as observed in the simulation. The 



phase shifts of those incidence soli tons are given by (3.7 1, 



AjCr 



1,4] 



-0.934, 



Axi 



[2,3] 



-1.518. 



Notice that the [1, 4]-soliton behind the [2, 3]-soliton shifts x^^ = x^^ ^-^ + Ajcji 4] = -0.307, and it gets a larger shift in 
the upper part. The large bending in the simulation is generated by a strong perturbation generated at the intersection 
point at f = 0. However, the bending part disperses and the intersection point gets away from this part. 



The bottom graph in Figure 5.12 shows the decay of the error. We emphasize that the error function gets a rapid 
change in the early times. This implies that a drastic change in the profile of the solution due to a large perturbation to 
the exact solution. However after f = 3 the solution has a monotone convergence to the exact solution, implying that 
the dispersive waves are getting away from the domain Dj. as we expected. 
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